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Abstract 

The exact calculation of network reliability in a probabilistic context has been a 
long-standing issue of practical importance, but a difficult one, even for planar 
graphs, with perfect nodes and with edges of identical reliability p. Many approaches 
(determination of bounds, sums of disjoint products algorithms, Monte Carlo evalu- 
ations, studies of the reliability polynomials, etc.) can only provide approximations 
when the network's size increases. 

We consider here a ladder graph of arbitrary size corresponding to real-life network 
configurations, and give the exact, analytical solutions for the all- and two-terminal 
reliabilities. These solutions use transfer matrices, in which individual reliabilities 
of edges and nodes are taken into account. The special case of identical edge and 
node reliabilities — p and p, respectively — is solved. We show that the zeros of 
the two-terminal reliability polynomial exhibit structures which differ substantially 
for seemingly similar networks, and we compare the sensitivity of various edges. We 
discuss how the present work may be further extended to lead to a catalog of exactly 
solvable networks in terms of reliability, which could be useful as elementary bricks 
for a new and improved set of bounds or benchmarks in the general case. 

Key words: network reliability, star-triangle transformation, transfer matrix, 
zeros of the reliability polynomial, sensitivity, algebraic structures 



1 Introduction 

Network reliability has long been a practical issue, and will remain so for 
years, since networks have now entered an era of Quality of Service (QoS). IP 
networks, mobile phone networks, transportation networks, electrical power 
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networks, etc., have become "commodities." Connection availability rates of 
99.999% are now an objective for telecommunication network operators, and 
premium services may only be deployed — and correspondingly billed — if 
the connection reliability is close enough to unity. Reliability is therefore a 
crucial parameter in the design and analysis of the various kinds of networks 
we daily use. 

Not surprisingly, the study of network reliability has led to a huge body of 
literature, starting with the work of Moore and Shannon [43], and including 
excellent textbooks and surveys [6,7,24,59,60]. In what follows, we consider a 
probabilistic approach, in which the network is represented by an undirected 
graph G = (V, E), where V is a set of nodes (also called vertices) and E is a set 
of undirected edges (or links), each of which having a probability p n or p e to 
operate correctly. Failures of the different constituents are assumed to occur 
at random, and to be statistically independent events. Among the different 
measures of reliability, one often considers the A;-terminal reliability, namely 
the probability that a given subset K oik nodes (K C E) are connected. The 
most common instances are the all-terminal reliability Kel^ (K = E) and the 
two-terminal reliability Rel 2 (s — > t), which deals with a particular connection 
between a source s and a destination t. Both of them are affine functions of 
each p n and p e . 

The sheer number of possible system states, namely 2' s ' + ' l/ ', clearly precludes 
the use of an "enumeration of states" strategy for realistic networks, and 
shows that the final expression may be extremely cumbersome. Consequently, 
most studies have considered graphs with perfect nodes (p n = 1) and edges 
of identical reliability p; radio broadcast networks have also been described 
by networks with perfectly reliable edges but imperfect nodes [1,32]. It was 
shown early on — see for instance the discussion in [24,73] — that the calcu- 
lation of /c-terminal reliability is #P-hard in the general case, even with the 
following simplifying and restricting assumptions that (i) the graph is planar 
(ii) all nodes are perfectly reliable (iii) all edges have the same reliability p. 
All reliabilities are then expressed as a polynomial in p, called the reliability 
polynomial. 

The difficulty of the problem has stimulated many approaches: partitioning 
techniques [26,78], sum of disjoint products [2,5,36,50,52,63], graph simplifica- 
tions (series-parallel reductions [43], delta-wye transformations [18,27,31,53,72], 
factoring [41]), determination of various lower and upper bounds to the reli- 
ability polynomial [6,8,13,15,20,24,57], Monte-Carlo simulations [29,39,44], ge- 
netic [21] and ordered binary decision diagram (OBDD) algorithms [42,51,76,77] 
The reliability polynomial has also been extensively studied [19,25,47], with 
the aim of deriving some useful and hopefully general information from the 
structure of its coefficients [19,25] or the location of its zeros in the complex 
plane [14]. 
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A recent breakthrough has been obtained by statistical physicists, who ob- 
served that the Tutte polynomial T(G, x, y) of a graph G is actually equiva- 
lent to the Potts model partition function of the g-state Potts model [61,74], 
which they were able to calculate for various recursive families of graphs [17]. 
The all-terminal reliability ReU(G,p) of such graphs is then deduced from 
T(G, 1, jz^)- Royle and Sokal [54] proved that the Brown- Colbourn conjec- 
ture [14] on the location of the zeros of the all-terminal reliability polynomial, 
while valid for series-parallel reducible graphs, does not hold for some fami- 
lies of graphs. While these results are extremely valuable to understand a few 
properties of graphs and all-terminal reliability polynomials, they still assume 
that nodes are perfect and that edges have the same reliability. 

In recent years, the tremendous growth of Internet traffic has called for a 
better evaluation of the reliability of connections in optical networks. This, 
of course, strongly depends on the connection under consideration. Actual 
failure rates and maintenance data show that a proper evaluation of two- 
terminal reliabilities must put node and edge equipments on an equal footing, 
i.e., both edge (fiber links, optical amplifiers) and node (optical cross- connects, 
routers) failures must be taken into account. The possibility of node failure has 
been considered in early papers [1,28,33], to quote but a few. Adaptation of 
algorithms to include imperfect nodes has been — sometimes controversially 
- addressed [40,45,70,71,77]. The two-variable approach for bounds to the 
reliability polynomial, by Bulka and Dugan [16] and Chen and He [20], is also 
worth mentioning. In order to be realistic, different edge reliabilities should be 
used too: for instance, the failure rate of optical fiber links is often assumed 
to increase with their length. 

In this work, we give the exact solution to the two-terminal reliability of a 
simple ladder network, displayed in Fig. 1, where successive nodes are labelled 
Si or Tj, and where the larger black dots mark the source s and terminal t we 
consider in our two-terminal reliability calculations (we choose these special 
nodes to lie at the graph extremities, since we can always reduce to this case 
by a series-parallel simplification). This network is a simplified description 
of a standard (nominal + backup paths) architecture, with additional con- 
nections between transit nodes for enabling the so-called "local protection" 
policy, which bypasses faulty intermediate nodes and/or edges. Such an ar- 
chitecture of "absolutely reliable nodes and unreliable edges," with up to 25 
edges, was chosen as Example 5 in [36] for a comparison of different "sum of 
disjoint products" minimizing algorithms, or by Rauzy [51] as well as Kuo and 
collaborators [42,76,77]. By letting the individual node and edge reliabilities 
take arbitrary values, we actually do not add to the complexity of the problem 
but make the internal structure of the problem more discernible (a similar 
approach has been fruitful in the context of graph coloring [10,12]). Indeed, 
it is then easier to exploit to the full the recursive nature of the ladder graph 
while using the delta- wye transformation for graphs with unreliable nodes [31]. 
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We show that the two-terminal reliability has a beautiful algebraic structure, 
as its exact expression is given by a product of 3 x 3 transfer matrices (see 
eqs. (13) and (16) below for the configuration of Fig. 1(a), or eqs. (17) and 
(18) for Fig. 1(b)). Consequently, it can also be determined for an arbitrary 
size (length) of the network. 



S G S 1 S 2 S n ., S n 




Fig. 1. Different source-terminal connections for lattice graphs. 

Our aim is (i) to give a detailed derivation of the final results, so that re- 
searchers or engineers involved in reliability studies can readily use an easy- 
to-implement formula, since worksheet applications are well up to the task 
(ii) evaluate the implications of these findings, which suggest a new and po- 
tentially fruitful approach to reliability in the case of more general graphs or 
in combinatorial problems such as the enumeration of self-avoiding walks in 
lattices of restricted width (iii) emphasize anew the importance of algebraic 
structures of the underlying graphs in the determination of network reliability 
[11,59]. 

Our paper is organized as follows: In Section 2, we briefly recall the formulae 
for the delta-wye transformation for unreliable nodes. In Section 3, we define 
the notations for the different edge and node reliabilities, used in the detailed 
derivation of the main results. We then consider in Section 4 the case where 
all edges and nodes have identical reliabilities p and p, respectively, the size of 
the network appearing simply as an integer n. We give there the analytic so- 
lution of the two-terminal reliability ReF^p, p;n) for different configurations, 
and their asymptotic form when n — > oo. In Section 5, we expand the preced- 
ing expressions for p and p close to unity, and determine when it is possible to 
approximate the two-terminal reliability in terms of multiple failures. We also 
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provide in Section 6 the associated generating functions, a very useful tool in 
combinatorics perfectly suited to reliability studies, since they encode all the 
necessary information in a beautifully simple form. Prompted by the nearly 
universal character of the Brown-Colbourn conjecture [14], we then address 
the location of zeros of the two-terminal reliability polynomials in Section 7, 
and show that their structures may well be quite distinct even for seemingly 
similar networks. We then give in Section 8 a short glimpse of the sensitivity 
issue [37,55], namely the influence of a given component to the overall reli- 
ability, by comparing the contributions of particular edges (the rungs of the 
ladder). For the sake of completeness, we derive in Section 9 the all-terminal 
reliability for the ladder networks under consideration, for arbitrary values of 
edge reliabilities. Finally, we conclude by indicating several directions in which 
the present results may be further extended, so that, for instance, a catalog 
of exactly solvable networks — in terms of reliability — may be given rapidly 
[66,67], which could be useful as elementary bricks for a new and improved 
set of bounds or benchmarks for alternative methods in the general case. 



2 Triangle-star transformation for unreliable nodes 



The triangle-star — also called delta-star, delta- wye, and A— Y — transforma- 
tion has been used many times to simplify calculations of network reliability 
[18,24,27,31,41,53,72], even though it has mainly been applied in a perfect 
nodes context, to provide upper and lower bounds to the exact reliability. 
Here, we fully exploit this transformation in the case of imperfect nodes in 
order to obtain exact results. Since it plays a crucial part in the derivation, 
we recall the formulae first derived by Gadani [31]. 

Consider three particular nodes A, B, and C of a network represented in the 
left part of Fig. 2 (the reliability of the nodes are given by the same uppercase 
variables, in order to avoid an unnecessary multiplication of notations). The 
reliability of the edge connecting A and B is given by (lowercase) c, with similar 
notation for the remaining edges of the triangle. The aim of the triangle-star 
transformation is to replace the triangle by a star, which is possible by the 
addition of a new unreliable vertex O and three new edges connecting O to 
A, B and C, with reliabilities pa, Pb, and pc, respectively. Both networks are 
equivalent provided that the following compatibility relations hold [31] 



PAOpc = b + acB — abcB, 
PaO pb = c + abC — abcC, 
Pb O pc = a + be A — abc A, 
PaO PbPc — a b + be + ac — 2 abc. 



(1) 
(2) 
(3) 
(4) 
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Fig. 2. Triangle-star transformation for unreliable nodes. A, B, and C are the node 
reliabilities, with a, b, and c the edge reliabilities of the initial network, and pa, Pb, 
pc, and O those of the transformed network. 

Note that the first three equalities correspond to the probability that the two 
nodes under consideration are connected, while the last gives the probability 
that all three nodes are connected. A word of caution — given by Gadani — is 
worth mentioning in the case of successive triangle-star transformations: the 
triangles should have no common edge or node. 



3 Derivation of the main results 

We first name the different edge and node reliabilities of the ladder diagrams 
represented in Fig. 1(a), and detail how we can use the triangle-star transfor- 
mation of the preceding section. Assuming that So is always the source node, 
and using node reliabilities Si or Tj, we note a n the reliability of the edge 
(S n -i, S n ), b n for (S n , T n ), and finally c n for (T n _i, T n ), as shown in Fig. 3(a). 



3. 1 Calculation of Rel 2 (So -> T n ) 

The two-terminal reliability can be written 



because the end-node reliabilities obviously factor themselves out of the ex- 
pression. The transformation from Fig. 3(a) to Fig. 3(b) is a mere series sim- 
plification, where a single edge of reliability a n b n S n is now connecting S n -i 
to T n . We can then apply the triangle-star transformation to the triangle con- 
stituted by {S^-i, T n _!, T n }, and define new edges pi, p 2 , p 3 , and a node O as 
represented in Fig. 3(c). From eqs. (l)-(4), we have 



Rel 2 (5'o -> T n ) = T n F n (6 , a 1 , b u c u ■ • • , a n , b n , c n ). 



(5) 
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Fig. 3. Two-step decomposition of the ladder network from (a) to (c), using the 
triangle-star transformation. 

p 1 Op 2 = a n b n S n + b 

n—l C71 Tfi— 1 & n bn—1 bn &n (6) 

PiOp 3 = c n + b 

n—i @-n bn S n —i S n a n b n —\ b n c n S n —i S n , (7) 

p 2 0p 3 = b 

n—l &n bn &n $n -^n bn—1 bn &n $n -^n> (8) 

Pi Op 2 ps = a n b n c n S n + 

&n bn—1 b n S n + b 

n—i £n 2 a n b n —\ b n c n 

Sn- (9) 

It is then obvious that the two-terminal reliability of the new network displayed 
in Fig. 3(c) must be formally similar to the one we started from. The product 
T n Pi O (a dangling link) is easily factored, and we are left with the calculation 
of Rel 2 (S'o — > O). This is nothing but the calculation of Rel 2 (So — > T n -i), 
where O has replaced T n -i and where adjustments are due: b n -i and c n -i 
must be replaced by p 2 and c n -iPzT n -i, respectively. In short, we have 

Fn = Pi O F n _i(6 n _i — > p 2 , ip 3 T n _i). (10) 
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This relation between F n and F n _i, in which the recursive nature of the ladder 
graph is now fully apparent, can be further simplified. Indeed, assuming for 
F ra the simple form 



f n @n bn S n ~\~ 'Yra C n -\- 5 n b n C n S n 



(11) 



in eqs. (5) and (10), and replacing all occurrences of p 1 , p 2 , ps, and O, we 
easily find 



\ 6n J 



M n ■ 



( 8 ^ 

Pn-l 
ln-1 
\ $11-1 J 



(12) 



with the transfer matrix M n given by 



M n = 







0"n bn—1 Cn— 1 -^n— 1 


b n -i Sn-i T n -i 




b n -l c n-l S n -l T n _\ 


—a n b n _i S n -i T n „i 


—a n b n _i c n _i S n -i T n _i 


a n (1 — 2 6„_i) C n _x S'n-i T n _x 



.(13) 



We now need to calculate Fi. This is straightforward, since we have a simple 
series-parallel graph: 



Fi = So (ai fei 5i + b ci T - 6 ai &i ci Si T ), 



(14) 



from which we deduce (5\ = a\ So, 71 = 60 So ^0, an d ^1 — —bo a i So T , or even 
more simply 



7i 
\ 5l J 



M 1 



(15) 



Note that eq. (15) is correct, even if c is undefined; c may be set to for the 
sake of simplicity. We finally obtain the beautifully simple and compact result 
for the two-terminal reliability between nodes S and T n : 



Rel 2 (So -> T n ; n > 1) = T n ■ (b n S n , c n , b n c n S n ) ■ M n M n -i ■ ■ ■ M x 



\ / 



■ (16) 
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with M n given by eq. (13). We can also write T n ■ (b n S n ,c n ,b n c n S n ) = 
(0,1,0) • M n+1 . Note that Rel 2 (S — > T ) = b S T Q may be obtained by 
replacing M n ---M 1 , which is undefined for n = 0, by the 3x3 identity 
matrix — or no matrix at all. It is also worth remarking that each edge or 
node reliability — which may be arbitrary — appears once and once only in the 
different matrices or vectors, respecting the general property that Rel 2 (So — > 
T n ) is an affine function of each component reliability. We have independently 
checked the above expressions with a sum of disjoint products method for 
the first values of n. Numerically speaking, the formula given in eq. (16) is 
very easy to implement, even in a worksheet application where multiplication 
of 3 x 3 matrices is routinely performed; it also easily applies for ladders of 
arbitrary size, which is also quite important. This is the final answer to the 
two-terminal reliability of Example 5 of [36], which only gave the number of 
disjoint products for up to 25 unreliable edges, and of benchmark networks 
#29-#30 of [42,51,76,77], where Rel 2 (S -> T n ) was calculated for n = 19 and 
n = 99 using OBDD algorithms. 



3.2 Calculation of Rel 2 (S ^ S n ) 



Let us now turn to the calculation of the two-terminal reliability between two 
nodes located on the same side of the ladder, as S and S n in Fig. 1(b). The 
recipe is the same as above: the edges b n and c n are merged with the node 
T n to give a new equivalent edge between T n -i and S n , of reliability b n c n T n . 
Here again, we can use the triangle-star transformation to derive new recursion 
relations, and transfer matrices appear. The calculations being straightforward 
and similar to those in the preceding subsection, we give only the final result 



Rel 2 (S -> S n ;n> 1) r„) • M„M„_! • ■ ■ Mi 



V / 



(17) 



with 





( c 


bn-1 Sn-1 T n -l 


dn-1 b n -lS n -l 


T ri -1 






O n -l b n -i C n S n - \ Tn-l 


C-n T n —1 


&n—l b n —l C n S n - 


i T n - 


i 




^ On— 1 b n — i C n S n —i T n _i 


~b n ^\ C n S n -1 Tn-l 


a n -l (1 - 2fe n -l) C n 


S n -i 


T n -i 



and the additional convention a — 1. Since S n • (a n , b n T n , a n b n T n ) = (1, 0, 0) • 
M n +i, eq. (17) also applies in the case n — 0, which gives Rel 2 (S'o — > So) = So, 
as it should. 
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4 Identical reliabilities p and p 



While eqs. (16) and (17) give the two-terminal reliability in the general case, 
we consider here the special case where edges and nodes have reliabilities p 
and p, respectively. Our exact results may prove useful because identical edge 
reliabilities are usually taken for granted, and combinatorial aspects play an 
important part in the corresponding calculations [24] . We show in this section 
that the two-terminal reliability of the ladder graph has a simple analytic 
expression for arbitrary n, which we derive from the above expansions. 

4.1 Calculation of Relf 0_>T " ) (p, p) 

Taking Oj = 6j = q = p, and Si = 7* = p, we get 





VP 


p 3 p 2 


\ n 

pV | 


(l 


Rei So ^\p,p) = (pp 2 ,pp,p 2 p 2 ) ■ 


VP 1 


VP 


p 2 p 2 


■ 




K -P 2 P 2 


_p3 p2 


p 2 (l-2p)p 2 ) 





We could of course replace (p p 2 ,p p,p 2 p 2 ) by (0, 1,0), provided that the matrix 
exponent is n + 1 instead of n. The next step is to find the eigenvalues and 
eigenvectors of the matrix appearing in eq. (19), so that the final result may 
be further simplified. Obviously, a p n p n factor can be factored out of the 
matrix contribution. We are left with the determination of the eigenvalues 
and possibly eigenvectors of 





( i 


p 2 p 


2 \ 
P p ' 




M T = 


P 


1 


VP 


(20) 




K -PP 


-p 2 p 


p{\ -2p)p) 





From the first two rows of Mr, 1 — p p is obviously an eigenvalue of the matrix. 
By using standard linear algebra, the characteristic polynomial of Mt can be 
factored, the three eigenvalues x and x± of M T being simply 



x = l-pp, (21) 

_ 1 + 2 p ( 1 - p) p ± V 1 + 4 p 2 p - 8 p 2 + 4 p 4 p 2 
x ±- 2 • { ) 

While the determination of the eigenvectors may be used to finish the calcu- 
lation, this is by no means necessary. Indeed, Mt being diagonalizable — the 
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three roots are distinct — Rel^^^Go, p) must be of the form 



From the first values of Relf °^ Tn) (p, p) for n = 0, 1, 2 



Ref ^ To) (p, p)=pp 2 , 
Relf°- Tl) (p,p)=p 2 p 3 (2-p 2 p 
Relf°^ T2) 



(p, p) =p 3 p 4 (3 - 2p 2 p + p 2 p 2 (1 - p) (1 - 2p)) , 



(24) 
(25) 

(26) 



which have been independently checked by a sum of disjoint products proce- 
dure, and the inversion of a 3 x 3 Vandermonde matrix based on xq and x±, 
we obtain a and a±. Finally 



Rel (5o-T„) 



(P, P) 



l-pp) n+1 + (l+pp) 



X 



n+l 



X 



n+l 



Xa 



X. 



-pp(l - 2p + pp) 



x + — x^ 



(27) 



Setting p = 0.9 and p = 0.9 or 1 for n — 19 or 99 confirms the numerical 
results found for the 2 x 20 and 2 x 100 ladder networks [42,51,76,77]. 

We discuss the implications of eq. (27) in section 4.4 below. We will see in 
Section 6 that the simplicity of the result stems from our choosing a unique 
reliability p for all edges. Before that, we turn to other two-terminal configu- 
rations. 



4.2 Calculation of Relf 0_>5n) (p, p) 



In the So — > S n configuration, eqs. (17) and (18) lead to 

/ 

Rei S ^ Sn \p,p)=p(p,pp,p 2 p) ■ 







1 1 


VP 


pp 2 


p 2 p 2 








P 3 P 2 


VP 


P 3 P 2 


— p 3 p 2 


—p 2 p 2 


p 2 (l -2p) p 2 



V 



1^ 



V / 



.(28) 



Note that while a® = 1 is the only Oj which should not be set to p, we can 
compensate this missing p term in Mi by introducing a 1/p term in the right- 
most vector in eq. (17). The procedure is the same as above: a (pp)™ prefactor 
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can safely be extracted from the matrix on the right-hand side of eq. (28), so 
that we have to find the eigenvalues of another matrix, namely M s 





( 1 


P 


VP 




P 2 P 


1 


v 2 P 




K -P 2 P 


-VP 


p(l-2p)p 



(29) 



Unsurprisingly, the eigenvalues of Ms are again given by xq and x±. Find- 
ing the prefactors associated with Xq and x± from the first three values of 

Re\2 So ~^ Sn \v, p) is straightforward, and the final expression is 



Rel (5 -5 n ) 



(P, P) 



pU pn+l 



-pp(l-2p + pp) 



X^ 



x_ 



(30) 



It is worth noting that the expressions given in eqs. (27) and (30) only differ 
by a '±' sign. 



4.3 Calculation of Rel^ ^' (p, p) 



Let us introduce another configuration — the symmetrical ladder network of 
Fig. 4 — which looks quite similar to the two previous ones because (i) it has 
been considered in the literature as a way to describe a "gamma multistage 
interconnection network" [50], an architecture demonstrating the potential 
usefulness of recursive algorithms by Hansler [34] or Fratta and Montanari 
[30] (ii) it will shed some light on what can be expected from the study of the 
two-terminal reliability polynomial, as regards the location of its zeros. 




Fig. 4. Symmetrical ladder network. 
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We do not give the exact expression of Relf °^ n> for arbitrary edge and node 
reliabilities, since it can be derived from eqs. (13) and (16) by merely taking 
To — S n — bo — b n — 1. We limit ourselves to the (p, p) configuration and find, 
after a treatment similar to that of the previous sections, 



n pn+1 
1+ — X- 



(2-p) (x n + -x n _) 
+p(l-2p + pp) (x™- 1 



x 



n-l 



(31) 



Note that x is absent here: only two eigenvalues appear, even though the 
network architecture differs from the previous ones by two edges only. 



4-4 Variation of the eigenvalues with p and p 



From eqs. (27), (30), and (31), it appears that the two-terminal reliabilities 
for the ladders are given by sums of terms such as Aq, A", and A", with 

\ =pp(l-pp), (32) 
^± = \pp (l + 2p(l -p)p± v /l + 4p 2 p-8p 3 p 2 + 4p 4 p 2 ), (33) 

where n appears only in the exponents. A natural first step is to assess these 
quantities as p and p vary, in order to assess the relative contributions of 
the different terms. For the sake of simplicity, we have displayed on Fig. 5 
the variation of Ao, A + , and A_ with p, when p — 1. All three quantities are 
positive and vanish when p goes to zero, as expected. However, only A + tends 
to 1 when p — * 1 (remember that for p — > 1, Rel2 — * 1, implying that one of 
the eigenvalues tends to 1). In this particular case, we have A + > A > A_ in 
the entire range < p < 1. 



An important issue is the behavior of ReV 2 n ' n ' n '(p,p) as the size of the 
network — or n — increases. We see on Fig. 5 that for p = 0.9, A + is about 
ten times larger than Ao and A_. Consequently, when n is sufficiently large - 
this may be valid for n ~ 10 already — we have to a very good approximation 

Rel { 2 So ^ Tn ' Sn ' Un \p, p) w a + A™ , (34) 



which is observed in the numerical evaluations of Kuo et al. and Rauzy 
[42,51,76,77]. Basically, adding one "cell" to the ladder network will merely re- 
sult in another A + factor in the two-terminal reliability, so that A + represents 
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0.2 0.4 0.6 o.a 1 

Fig. 5. Dependence of Ao, A + , and A_ on p, when p = 1. 

a scaling of the reliability; it is as if we had n equipments of reliability A + in 
series. When p ^ 1, this general behavior is the same, even if A_ is negative for 
2^— < V < 1, A + still prevails, and the power-law behavior of eq. (34) holds. 



5 Failures decomposition 

In some instances, when p and p are assumed to be very close to 1, calculations 
often involve the unreliability of the connection [/ 2 = 1 — Rel 2 . Let q = 1 — p 
and i] = 1 — p. The network unreliability is then expanded in q and rj, even 
though one should be careful about such a practice, when q and i] differ by 
orders of magnitude. For instance, keeping only three terms in [/ 2 , we find 
from eq. (27) that the following formula is valid for n greater than 3 (but not 
too large, of course) 

U 2 {S ^ T n ) = 2r] + {3n - A) r] 2 - (12n - 16) r] 3 + 4npq 

-(20 ra - 6)r] 2 q + (n + 2) q 2 - (4n + 14)r/g 2 

+2 (n - 1) q 3 + • • • (35) 
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We must keep in mind, however, that q and r\ must be sufficiently small for 
such an expansion to make sense. Limiting ourselves to rj = (the perfect 
nodes configuration), we get 

n 2 4- 7 n 4- 1 f) 

U 2 (S ^T n ) = (n + 2) q 2 + 2(n-l) q 3 - n + ^ + g 4 + ■ ■ ■ (36) 

with the same result for U 2 (S — > whereas 

f/ 2 (5o-^)=ng 2 + 2(n-l)g 3 -^±^)g 4 + --- (37) 

The "small parameter" in the expansion is therefore of the order of nq. A 
rule of thumb is to ensure this quantity is small indeed when considering 
expansions such as eq. (35). 



6 Generating functions 



There is yet another way to express all the abovementioned results is a simple 
and compact form, namely that of the so-called "ordinary generating func- 
tions", which are a fundamental tool in combinatorics [65] and extensively 
used in [17]. The matrix transfer formalism and results obtained in the pre- 
ceding sections clearly indicate the existence of recursion relations between 
successive reliability polynomials (see [10] for an interesting counterpoint for 
the chromatic polynomials of graph ladders), which we can write in an even 
more compact form than eqs. (27) and (30). If we can write — assuming all 
the eigenvalues are distinct — 

Rel ( 2 n) (p,p)= £ a,i A™, (38) 

ie{o,+ ,-} 

the ordinary generating function Q(x) defined by 

oo 

g(x) = Y,^t\p,p)x n (39) 

n=0 



is easily calculated, since we have to sum geometric series : 



3 



Af(x) 



' w -J£_,^-i&- (40) 
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Of course, the summation can be performed too when oli is a polynomial func- 
tion of n. Q (x) is a rational fraction of x (see [65] for a thorough discussion of 
the topics). Its denominator T>{x) = _ \ x) m ay be deduced very simply 
from the characteristic polynomial P ca .ra,c{. x ) of the unfactored transfer matrix 
by V(x) = (— x) 3 P ca , Tac (l / x) . The determination of G(x), which encodes all 
the two-terminal reliabilities for the recursive family of graphs indexed by n, 
is actually quite simple: we only need to multiply the known leading terms of 
the expansion of Q(x) by V(x) to see J\f(x) emerge. 

For instance, in the case of Rel 2 S ' ~ + ' rn ' ) (p, p), we easily find V T (x) from eq. (19) 



V T (x) = l-ppx (2 + pp (1 - 2p)) + p 2 p 2 x 2 (1 +pp(l - 2p) (2 -pp)) 
-p A p 4 x 3 (1 - pp) (1 - 2p + pp) , (41) 

and deduce from eqs. (24)-(26) the simple result (mathematical softwares 
allowing symbolic calculus make such computations an easy task) 

M T {x) =pp 2 (l - (1 -p)p 2 p 2 x) . (42) 



The preceding treatment works well too if all the aj, 6j, and q are replaced by 
a, b, and c, respectively, instead of by a unique value p. The final result reads, 
after simplification 



2 N T (a,b,c,p;x) 



G T (a,b,c,p;x)=bp ' ' ' ' (43) 

V T (a,b, c,p;x) 

A/r(a, b, c, p; x) = 1 — a (1 — b) cp 2 x, (44) 

T> T (a, b, c, p; x) — 1 — [a + c + acp — 2abcp] px 

+a c [1 + (1 - 2 6) (a + c) p - (1 - a - c) b 2 p 2 ] p 2 x 2 
-{l-bp)(l-2b + bp)a 2 c 2 p i x 3 . (45) 

When a, 6 and c are not identical, no simple factorization of T>t occurs, so 
that we must use the roots of a polynomial of degree three in x, the literal 
expressions of which are not concise. However, the expression of the generating 
function makes it possible to calculate the two-terminal reliability by a mere 
partial fraction decomposition, which can be handled quite routinely. 

In the case of Relf°^ 5n) , Q s {x) = M s (x)/V s (x), with V s (x) = V T (x) and 



Af s (x)=p-pp 2 (l+pp(l-2p)-p 2 p 2 (l-p)) x 

+P 3 P 4 (l-2p + p 2 p(2-p)) x 2 . (46) 



16 



While this expression seems more complicated than eq. (42), we can check 
that 

Qs(x) - Gt(x) = - j- r— , 47 



which could be directly deduced from eqs. (27) and (30). 
The generating function Q v for S — > U n is given by 

gu (x) = pp 2 x (2-p)+p 2 p(l-2p + pp)x 

yuK> Pi l-pp(l + 2p(l-p)p)x + p :i pZ(l-2p + pp)x 2 ' { 1 



Here, the denominator is only a polynomial in x of degree two, leading to the 
decomposition in eq. (31). We have included the case n — 1, for which the 
reliability is that of two links in parallel between So and Lq, namely p 2 (2p — 
p 2 ); this has the advantage of simplifying the rational Gu{x), and making 
transparent the derivation of eq. (31). 



7 Zeros of the two-terminal reliability polynomials 



One way to understand the structure of the different reliability polynomials is 
to study the locations of their zeros in the complex plane. Such a study has 
been fruitfully performed in the case of the chromatic polynomial [10,12,56], 
most notably in the context of the four-color theorem. In the reliability con- 
text, some effort has been done to discover general properties for the all- 
terminal reliability ReU(p) [19,25,47], its main byproduct being the Brown- 
Colbourn conjecture [14], according to which all the zeros are to be found 
in the region |1 — p\ < 1. Although valid for series-parallel graphs, this re- 
markable conjecture does not strictly hold in the general case (but not by far) 
[54]. As mentioned in the introduction, the all-reliability polynomial is linked 
to the Tutte polynomial, an invariant of the graph. It has also been studied 
extensively by Chang and Shrock for various recursive families of graphs [17], 
who give the limiting curves where all zeros of the polynomials converge. 

In this section, we want to locate the zeros of the two-terminal reliability 
polynomial, in order to see whether some insight may also be found in this 
case. Truly, this polynomial depends on the couple (source, terminal), but 
some structures could still be expected. We show indeed that the zeros tend 
to aggregate along sets of curves, which can substantially differ even for two 
seemingly similar ladder graphs. 
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7.1 Calculation of the limiting curves: General results 



As n grows, the number of zeros of the reliability polynomial in the complex 
plane increases. Because of the matrix transfer property, we have recursion re- 
lations between polynomials corresponding to successive values of n, and the 
solution looks like eq. (23). The general treatment of the problem has been 
done by Beraha, Kahane, and Weiss [9], but may be understood in the fol- 
lowing, simplifying way: if the reliability polynomial has the form J2i a i \ {p) n 
(where \ are the eigenvalues of the transfer matrix), then at large n, only the 
two eigenvalues of greater modules, say Ai and A2, will prevail, so that the 
reliability polynomial will vanish when |Ai(p)| = | A2 Cp) | (of course, it might 
be three or more eigenvalues of equal moduli; our simplification is sufficient 
here) . This defines a set of curves in the complex plane, where all zeros should 
accumulate in the n — > 00 limit. The interested reader should refer to the 
work of Salas and Sokal [56] for a very detailed discussion of the convergence 
to the limiting curves. This behavior would not be modified when one or more 
of the ctj's is a polynomial in n [9]. 

In the following we limit ourselves to the perfect-nodes case (p = 1) and look 
for the limiting curves for the architectures considered in Section 4, because it 
is illustrative enough, and because all expressions are then simple and analyt- 
ical. The trivial solution p = will not always be considered in the following. 

7.2 Calculation for the symmetrical ladder (So — > U n ) 

Let us first consider the simplest case for Rel 2 («So — > U n ) given in eq. (31) 
because only two eigenvalues are present (see Fig. 4). We have to find p such 
that \x+/x-\ = 1, or 



1 



l + 2p(l-p) + v /l + 4p 2 (l-p) 2 



(49) 




Introducing t 



2p(l 



p), we must find t such that 



i + t + VT+¥ 



9 



(-7T < e < 




(50) 



1 +t- vm^ 



which can be rewritten as 




(51) 
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Squaring both sides of eq. (51) implies 
t 2 + 2 sin 2 ^ + 1 = 0, 



(52) 



which yields two possible solutions, t±, of modulus 1: 



6 / 6 

t± = -sin 2 -± iWl-sin 4 -. (53) 
2 V 2 



Only the solutions that verify eq. (51) must be kept ; they may be rewritten 
in condensed form 

«=-«* {-i (54) 



which represents a half-circle (the real part of t± must be negative), so that 
the searched values of p are finally 



P= 2 (55) 



It is then straightforward to give the coordinates of some particular points 
of the curves. On the real axis, the solutions for ip = are (1 ± \/3)/2, 
i.e., +1.366 and -0.366. The endpoints of the two curves are actually \ (1 ± 

yj^fi) ±i^^T, or 1.136 ± 0.393 i and -0.136 ± 0.393 i. 
For n = 10, we have 



Rel 2 (5 -> C/io) =P 10 (2 + 18 p + 68 + 100 p 3 - 134 p 4 - 746 p 5 - 648 p 6 
+ 1824 p 7 + 3818p 8 - 2354 p 9 - 10861p 10 + 2586p n 
+23080 p 12 - 7904 p 13 - 48624 p 14 + 79008 p 15 
-58432 p 16 + 24064 p 17 - 5376 p 18 + 512 p 19 ) . (56) 



We see on Fig. 6 that for n = 10 already, the zeros are located very close to 
the two well-separated arcs of quasi-circular shape defined by eq. (55). Note 
that with respect to the diagrams by Chang and Shrock for the all-terminal 
reliability polynomials calculated for several recursive families of graphs [17], 
we now have accumulation points in the region Re(p) < 0. When p ^ 1, the 
global appearance is unchanged, even if the zeros move away from the origin. 
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1.5 



1 




-0.5 



Fig. 6. Location of the zeros of Rel^^^ip): the dots corresponds to n = 10 and 
eq. (56), the continuous curves to the limit n — > oo given by eq. (55). The trivial 
solution p = is not displayed. 



When p approaches zero, the structure distorts itself to look like a circle of 
radius of the order of (2 p) -1 ^ 2 centered at (|, 0). The convergence of the roots 
to the limiting curve follows the general expression of Salas and Sokal [56] . 



7.3 Calculation for Rel 2 (S , — > S n ) and Rel 2 (5 , — > T n ) 



In these configurations, we now have three eigenvalues. We must therefore 
consider |x±/xo| in addition to \x + /x-\, calculated in the preceding section. 
We must find all the p's such that 



^H^jHW^,, (_<,<,), (57) 
Xq 2(1— p) 



so that 

Z = ±y/l + 4p 2 (1 - pf = 2 (1 - p) e ie - (1 + 2p (1 - p)) . (58) 
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Here again, we square both sides of eq. (58), keeping in mind we will have to 
check that the real part of the right- hand- side of eq. (58) must comply with 
our determination of the square root (our angle argument being taken between 
—n and +71, we choose Re^/ 7 ^ > 0). After simplification, we get 



l + 2p(l-p) = (1-p) e ie +pe~ 



(59) 



This equation has two possible solutions, namely 
p± = i (l + i sin 6 ± ^i + tl-costf) 2 " 



(60) 



Which values of 9 satisfy with eq. (58) ? Using eqs. (58) and (59), we find 

Z = (1 -p±)e ie -p± e~ ie = =Fcos0^/l + (l-cos0) 2 + i sin9 (1 - cos 9), (61) 



so that Re(Z) and Re(=)=cos6 l ) have the same sign. Since for x + /xo, we must 
have Re(Z) > and thus Re(=Fcos0) > 0, we find 



p+ 

< 


(f 


< 


1^1 


< 7T) 


V- 

V 


(0 


< 




<f) 


f 

V- 

< 


(I 


< 




< 7T) 




(0 


< 


1^1 


<f) 



(62) 
(63) 



We also have to check, when considering p[|x + /xo| = 1] for instance, that the 
third eigenvalue x_ has a smaller modulus. Consequently, using eqs. (59) and 
(61), 



'1-p) e id +pe~ id + ((1-p) e ie -pe 



[l-p] e ie +pe~ ie - {(1-p) e ie -pe- ie ) 



1 — p 



P 



>1, (64) 



and therefore that Re(p) < \, where p stands for p±. We deduce from eq. (61) 
that the p + solution is to be discarded. The same procedure for p[|x_/x | = 1] 
shows that only p_ should be considered. 

Finally, we must assess the modulus of x Q along the curves given by eq. (55). 
Here, the situation is slightly more complicated, and we checked numerically 
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that the p_ solution of eq. (55) should be discarded. The final result 
1 + Vl + 2e^ 



P 



\ (l - ^l + fl-costf) 2 + i sin 6) 



(-1 <^<f) 

(-7T < < 7T) 



(65) 



is displayed in Fig. 7. 



1.5 



1 




\ 




-0.5 



- — i — ■ — ' — ( — ■ — i — i — ■ — ■ — i — = — i — i — ' — ■ — >- 



-0,5 



0.5 



1.5 



Fig. 7. Location of the zeros of Rel^ (p): the dots corresponds to n = 10 
(eq. (66)), the continuous curves to the limit n — > oo given in eq. (65). 



When n = 10, we have 



Rel 2 (So -> Tio) =p 11 (11 + 155 p 2 -99p d + 40 p 4 - 907 p b - 296/ 
+ 1448/ + 3121/ - 1102 p 9 - 7989 p 10 - 1747p n 
+ 14806/ 2 + 4776p 13 - 24168 p 14 + 176p 15 + 35072/ 6 
-38016/ 7 + 19072/ 8 - 4864 p 19 + 512 p 20 ). (66) 

We have plotted in Fig. 7 the zeros for n — 10 and the n — > oo limit. Note that 
the asymptotic limit seems to be reached faster in the Re(p) > half-plane. 
The new curve endpoints of Fig. 7 are \ (1 - V2 ± i) sa -0.207 ± 0.5 i, and 
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I (1 — ~ —0.618. For p 7^ 1, the structure of the zeros spreads out from 
the origin, as in the preceding subsection. 

We want to emphasize that while the structure of the zeros in Figs. 6 and 
7 is similar in the Re(p) > half-plane — the n — > oo limit is identical - 
this is certainly not the case in the Re(p) < half-plane, even though the 
two graphs differ by two links only. Inferences from these structures of zeros 
should therefore be considered cautiously. 



8 Sensitivity 



An important parameter appearing in reliability studies is the sensitivity Sj 
[55] , which measures the influence of a particular equipment i (node or edge) on 
the global reliability; the sensitivity analysis is akin to the various importance 
criteria compiled by Henley et Kumamoto [37]. It is usually defined as the 
derivative of the total - - be it two- or all-terminal — reliability with the 
individual reliability pi 

QD 1 

s t = = Rel 2 ,Afe = 1) - Rel 2A ( Pt = 0). (67) 



The two definitions are equivalent, because Re^A is an affine function of all 
the pj's. The calculation can be performed quite easily, since the desired Pi 
appears in only one of the transfer matrices of Section 3. For instance, let us 
consider the sensitivity of the successive "rungs" of the ladder — a similar 
work could be done for the nodes, without difficulty In the following, we 
restrict ourselves to the case Rel 2 (S , — > T n ), with perfect nodes and identical 
edge reliabilities p, and want to assess the influence of the edge location to the 
overall performance of the network. 

From eq. (13), taking the derivative of a matrix with respect to a given b no 
is not difficult, since each of its elements is an affine function of b m . A new 
matrix appears in the problem, and when all edge reliabilities are set to p (we 
take p = 1 for the sake of simplicity), we have 
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p 3 


p 3 ) 


n—no 
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p 2 ^ 


^ = (0,1,0) • 
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K -p 2 
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V 


p3 


p 3 \ 


no 


M 


V 


V 









—p 2 


_p3 


p 2 (l-2p) j 







(68) 



For symmetry reasons, the contributions of fej and 6 6 _j are identical. For n — 6, 



s bo = (1 - p)p 6 (1 + p + 15 p 2 + 4p 3 - 18p 4 - 55p 5 + p 6 + 116p 7 





+24p 8 - 200p 9 + 144 p 10 - 32p n ), 








(69) 




= (1 -p) 2 p 6 (l + 2p + 16p 2 + 15 p 3 - 
+92 p 7 + 40 p 8 - 96 p 9 + 32 p 10 ), 


lip 4 


-60p 5 


-28p 6 


(70) 




= (1 -p) 2 p 6 (l + 2p+ 16p 2 + 14 p 3 - 
+88 p 7 + 40 p 8 - 96 p 9 + 32 p 10 ), 


14 p 4 


-61p 5 


- 20 p 6 


(71) 


Sb 3 


= (1 -p) 2 p 6 (l + 2p+ 16p 2 + 14 p 3 - 
+88 p 7 + 40 p 8 - 96 p 9 + 32 p 10 ). 


15 p 4 


-60p 5 


-20p 6 


(72) 
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P 

Fig. 8. Variation of the sensitivity Sens2 of the ladder's rungs, depending on their 
location, in the Rel2(So — ► Tq) case. 



Figure 8 shows that these sensitivities vary with p and that they uniformly 
decrease from the ladder's ends to its center. Note that only a few coefficients 
of the polynomials differ. If the length of the ladder increases, the maximum 
sensitivity will actually decrease in magnitude, and be located at higher values 
of p. Besides, if n is odd, then the minima are shared by the two central rungs, 
i.e., n = (n ± l)/2. 
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Because of the form of eq. (68), the final result is expected to be a sum of 
different power-law contributions, with products of two of the different eigen- 
values already observed for the reliability. The easier way is then to calculate 
the generating functions in a So — > T 2n ladder, for both b and b n = 6 cen trai- 
We finally get Q sens (bi) = E n =o s bi (S -> T 2n ) x n with £ scns = A/" S en S /^>sen S : 



Mens (Central) = 1 ~ p 2 (1 + p + 3 p 2 - 6 p 3 + 2 p 4 ) X 

+ (l-p) p b (1 + 3p - 3p 2 - 2p 3 + 2p 4 ) x 2 
V(l-P) 3 * 3 , (73) 

£>scns(&central) = (l ~ p 2 (1 ~ pf x) (l - p 3 (1 - P ) x) 



x (l - p 2 (1 + 2p + 2p 2 - 8p 3 + 4p 4 ) x + p 6 (1 - pf x 2 ). (74) 

whereas 

n s 1 - p 2 (1 + 3p 2 - 5p 3 + 2p A ) x + (2 - p)p & (1 - pf x 2 

^sens(»o) = 7 

(1 -p 2 (1 -pfx) (1 -p 2 (1 + 2p + 2p 2 - 8p 3 + 4p 4 ) x + p 6 (1 -pfx 2 ) 



The partial fraction decompositions of £/ S ens(fro) and ^ sen s (Central) lead to 



2 [ 2 ^ ^1 + 4^(1 -p) 1 

x 2n ^ 1 - 2p 2 

1 + 4p 2 (l - pf 



(76) 



and 



^ ^central 



i 2n \m „\2n , 4p 2 (1 - p) 

2 P 



U\-pf n + 








+ 4 + i 


x 2n | 






+ 4 1 



l + 4p 2 (l- 




l-2p 2 




^l + 4p 2 (l- 


-pf J 


1 — 2p 2 


\ 



so that the ratio Sf, contral /s& , m the large n limit, is given by 



2 



Sfeo 2 ^ v /l + 4p 2 (l-p) 2 / ' 



(77) 



s 6 ccntral 1 / x + 1 2p ^ (7g) 
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as displayed in Fig. 9 (for p = 1 — q close to 1, this ratio expands as 2q — 
6q 3 + 8g 4 + • • •)• 




Fig. 9. Ratio s bccntral /s bo for Rel 2 (S , -> T 2n ) in the large n limit. 



Actually, this limit is reached very quickly, since for the So — > T 6 ladder, the 
ratio Sb 3 /sb is barely distinguishable from the asymptotic limit in the entire 
range < p < 1, which is displayed in Fig. 9. It also shows that the overall 
reliability is more sensitive to "access" rungs located at the edges of the ladder 
network than to those in the middle, even if all edge reliabilities are identical. 



9 All-terminal reliability 



As mentioned in the introduction, the all-terminal reliability is another useful 
measure of the network availability, namely the probability that all nodes are 
connected. We refer the reader to the work of Chang and Shrock [17] for the 
explicit expressions of ReU for various recursive families of graphs, among 
which our simple ladder graph, when all edges have the same reliability p. For 
the sake of completeness, we show that their result may be slightly generalized 
for edges with distinct reliabilities. Here again, the final, analytical expression 
can be written in a concise form using transfer matrices. In the context of 
graph theory, this can be viewed as the factorization of a particular value of 
the multi-variate Tutte polynomial considered by Wu [75] and Sokal [64]. 
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Fig. 10. Successive deleting/contracting operations lead to a recursion relation be- 
tween lZ n and lZ n -i (see text). 

We expect the all-terminal reliability TZ n for the simple ladder diagram of 
Fig. 1 should exhibit the same behavior, with a transfer matrix depending on 
a n , b n , and c n , and a generating function that is a rational fraction of p (the 
common reliability of links) and x. Actually, the calculation of the all-terminal 
reliability may be somewhat easier because all nodes can be considered perfect 
without loss of generality [24]. We can then use the usual factoring decom- 
position to establish a relationship between lZ n and 7Z n -i- From Fig. 10, we 
easily find that 

K n = [(1 - a n )b n + a n (l - &„)] c n Tl n -\ + a n b n T^n-i{b n -i -> b n -i / / Cn). (79) 
Using the ansatz 1Z n = a n (a n + b n ) + f3 n a n b n , we get 



' (an-! + b n -i) C n 


^n— 1 b n ~\ C n \ 




K (a n -i + b n -i) (1-2 Cn) + c n (1 - 2 6 n _i) 


a n -i (c„ + 6„_i (1 - 3c n )) y 


y 



,(80) 



where actually ao = 1 and /3q = 0. Calling M n the matrix appearing in eq. (80), 
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we obtain (with the condition ao = in Mi) 
K n = (a n + b n , a n b n ) ■ M n ■ M n _ x ■ ■ ■ M 1 







When all the reliabilities are equal to p, the transfer matrix is simply 

\ 



M. 



R 



2p 2 


p3 


p(3- hp) 


p 2 (2-3p) i 



(82) 



The eigenvalues of Mr, £± = p 2 
tained, and finally 



4 - 3p± V12-20p + 9p 2 



Tlr, 



pV 12 -20p + 9p 



(/-n+l _ fn+l\ 



, are easily ob- 



(83) 



and 



G R (x) 



p 



1 — p 2 (4 — 3p) x + p 4 (1 — p) x 2 



(84) 



The last results were obtained in [17]: the polynomials aggregate at the curve 
defined by p = 1 + \e lip 1 with cosy? > |. It is worth stressing that the new 
scale C+ for the all-terminal reliability is different from the eigenvalues found 
in the two-terminal study, but that a power-law behavior is still the rule. 



10 Conclusion and perspectives 



We have found the solution of the two- and all-terminal reliabilities for a simple 
ladder graph, which corresponds nonetheless to realistic network architectures, 
especially in telecommunication networks for IP transport, but not only. Node 
and edge failures are put on an equal footing, and the simple formulae rely- 
ing on transfer matrices may be directly implemented, even in a worksheet 
application. When identical reliabilities p and p are assumed for edges and 
nodes, respectively, we have given the analytical solution of the two-terminal 
reliability, for which only numerical determinations were previously available. 
We have shown that, while it may be useful to study the location of the ze- 
ros of the two-terminal reliability polynomial, it may strongly depend on the 
existence of a limited number of particular edges, especially in the Re(p) < 
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region. We have also given the generating functions of the reliability polyno- 
mials, and provided a glimpse of sensitivity studies, which could be performed 
straightforwardly. 

Although we used a delta-star transformation to solve a particular family of 
graphs, the present work may clearly be extended in several directions, which 
we outline in the following, and which we shall develop elsewhere. 

10.1 Brecht-Colbourn ladder 

Another simple ladder graph is the Brecht-Colbourn ladder [13,48], which has 
been considered study, both for the two- and all-terminal reliability, 

in order to evaluate the quality of bounds for the reliability polynomial in 
the case of perfect nodes. It has been also brought forward in a special case 
of broadcasting network [32]. Using the same delta-star transformation, the 
application of which is slightly more tricky, we have been able to solve this 
problem exactly, where the intermediate nodes have a connectivity degree of 4 
instead of 3. The complete results — the transfer matrix is then of dimension 
four for imperfect nodes, but still of dimension three for perfect nodes; a few 
eigenvalues may be complex — will be given elsewhere [66], along with results 
for a generalized fan [3,46]. 

10.2 More general recursive families of graphs 

Even though the delta-star transformation may not apply successfully to all 
networks, it seems quite clear that a similar decomposition through transfer 
matrices of two- and all-terminal reliabilities should exist too for a ladder of 
K 4 graphs, or ladders of greater width. In order to make such calculations 
useful for applications, imperfect nodes as well as imperfect edges should be 
considered. All is needed is a recursion relation between successive graphs, 
when one "elementary cell" is added. This implies a new expression for the 
deletion-contraction theorem, in which the — most useful — linearity with 
respect to all individual edge or node reliabilities must be preserved, whereas 
to our knowledge the edge reliabilities have been renormalized to account for 
the unreliability of the nodes they connect [70,71]. This new expression will 
be given elsewhere [67], along with an application to other recursive families 
of graphs (K 4 ladders and K 3 cylinders). 

What should we expect ? Basically the same kind of behavior as detailed in 
the present work, with a factorization of the reliability in terms of transfer 
matrices, the dimension of which may substantially increase to reflect the in- 
terplay of different edges/nodes in the overall reliability, and the underlying 
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algebraic structure of the graph. It should be clear, however, that all graphs 
might not behave that way [10]. For instance, the generating function of the 
two-terminal reliability for the complete graphs K n , for which we give a recur- 
sion relation in the case of imperfect nodes and edges in the appendix, might 
not be a rational fraction. 

10.3 Failure frequency of systems 

Another performance index for networks, mostly investigated in the context of 
power distribution systems [35,58,62], is the failure frequency of a connection. 
Calculations can readily be performed when the assumption is made that all 
equipments fail at a constant failure rate A e and brought back to operation 
via a constant repair rate /i e (and p e = — ^ e ). The failure frequency is then 

the sum over all equipments of the products of the individual failure rate, 
reliability and sensitivity (as defined in Section 8). The present results clearly 
indicate that such calculations are made very simple by the transfer matrix 
formulation, even for extended networks. We shall present various examples 
in a forthcoming paper [68]. 

10.4 Reliability incertitude 

In recent years, Coit and collaborators [22,23,38] have considered the possible 
influence of the reliability uncertainty to the determination of an "optimal 
network architecture" , and have mostly dealt with series-parallel architectures. 
In the context of meshed networks, and recalling our assumption of statistically 
independent failures for the network constituents, we could safely replace in 
the final expression of the reliability each p e or p n by the random variable 
P e or P n of which it is the mean value. We can thus straightforwardly assess 
the influence of the full statistical properties of P e or P n (variance, skewness, 
kurtosis, etc.). This issue will be reported elsewhere. We can also replace, for 
non-reparable systems, each p e by the often used exp(— \ e t). 

10.5 Bounds, combinatorics, etc. 

There are obviously many directions in which this work may be useful in the 
determination of bounds. The first one is to use our present results on simple 
ladder graphs, or their future extensions [66,67], as possible upper or lower 
bounds to more complex graphs. If a particular graph looks like a special 
instance of a recursive family of graphs, we expect the generating function 
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to be a rational fraction again. The dimension of the corresponding transfer 
matrix may be probed by trying to find recursion relations between succes- 
sive reliability polynomials. Of course, if the dimension of the corresponding 
transfer matrix is large, the degree of the numerator and denominator of the 
fraction may be too large for a complete solution to be obtained easily. Even 
so, the knowledge that the generating function is rational may be an indica- 
tion that an approximate generating function might be quite useful, since Pade 
approximants [4] are known for their devilish knack of getting very close to 
the exact function [49]. This suggests that we may profitably consider rational 
fractions — deduced from approximate generating functions — for bounds, 
instead of exclusively relying on polynomials. This will be demonstrated in 
detail elsewhere [69]. 

Our calculations may also provide some information on some combinatorial 
issues, such as the enumeration of self-avoiding walks on lattices of restricted 
width. Finally, the exact results found for classes of arbitrarily large net- 
works may prove useful for testing different algorithms (Monte Carlo, genetic, 
OBDD, etc.) in numerically exacting configurations, where edge and node 
unreliabilities have to be taken into account. 
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A Two-terminal reliability for a complete graph with imperfect 
nodes 



The exact all-terminal reliability polynomial A n = ReU(-Kn) for the complete 
graph K n with perfect nodes and identical edges has been known for decades, 
as well as the corresponding two-terminal reliability polynomial T n = Ke\2(K n ) 
([24], pp. 33-34). They may be found by recurrence, starting with A± — 1: 



4. = 1 " E 



n-1 ( n -\\ 



A^l-p) 



j (n-j) 



(A.1) 



and 



E 

i=2 



n 



Aj(l-p) j ^- j \ 



(A.2) 



Taking imperfect nodes with identical reliability p into account does not sub- 
stantially change the all-terminal reliability, since one merely has to multiply 
the result of eq. (A.l) by p to the power of the number of operating nodes [24]. 
The situation is quite different for ^-terminal reliability, and in particular for 
two-terminal reliability, since all the possible paths between the source and 
the destination may visit different numbers of nodes. This can be very useful 
since the complete graph configuration represents the highest possible upper 
bound, because all its nodes are connected. 

Defining T n (p, p) the two-terminal reliability polynomial with imperfect edges 
and nodes — so as to keep the notation of Colbourn [24], we have 



T 2 (p, p)=PP 2 , 

T3(p,p)=pp 2 +p 2 p 3 -p 3 p\ (A.3) 
T*(p, p)=pp 2 + 2p 2 p 3 - 7p 4 p 4 + 7p 5 p 4 - 2p 6 p 4 + p 3 (-2 p 3 + 2 p 4 ) , 

and more generally 

P j=l \ j — 1 I 



with Vi(p) = p, V 2 (p) = P 2 -p 3 , and V 3 (p) = 2p 3 - 7p 4 + 7p 5 - 2p 6 , etc. We 
can then invert eq. (A.4) given for T„(p, p = 1) to obtain all the Pj(p)'s, using 
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the inversion formula 



P;(P) = E (-1) 



fc=i 




T fc+1 (p,p = l). 



(A.5) 



Combining Eqs. (A. 4) and (A.5) gives an easily computable expression for 
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